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Abstract. We study the dynamics of a kink in a one-lane asymmetric simple 
exclusion process with detachment and attachment of the particle at arbitrary sites. 
For a system with one site of detachment and attachment we find that the kink is 
trapped by the site, and the probability distribution of the kink position is described 
by the overdumped Fokker-Planck equation with a V-shaped potential. Our results can 
be applied to the motion of a kink in arbitrary number of sites where detachment and 
attachment take place. When detachment and attachment take place at every site, we 
confirm that the kink motion obeys the diffusion in a harmonic potential. We compare 
our results with the Monte Carlo simulation, and check the quantitative validity of our 
theoretical prediction of the diffusion constant and the potential form. 



1. Introduction 

An asymmetric simple exclusion process (ASEP) is a simple nonequilibrium stochastic 
lattice model of transport process. ASEP is defined on a one- dimensional lattice, and 
particles hop to the neighbouring site when the site is empty. The transportation of 
particles is induced by the asymmetric hopping of the particles. Recently much attention 
is paid to ASEP [1, 2, 3], not only because there exists the exact solution under the 
open boundary condition [4, 5, 6, 7], but also it is applicable to various transportation 
phenomena. In particular, the uni-directional ASEP which is called totally asymmetric 
simple exclusion process (TASEP) has been studied extensively, because (i) TASEP 
is the simplest ASEP and (ii) TASEP keeps the essence of nonequilibrium transport 
processes such as the exclusion interaction between particles and the drift of particles. 
In fact, TASEP may be regarded as a simplified model of traffic flow [8]. 

We can also apply ASEP to biological problems. TASEP is first introduced as a 
model to explain the process of creation of the messenger RNA [9]. The model with 
detachment and attachment on each site in TASEP with open boundary conditions is 
known as PFF model [10], named after the authors' names: Parmeggiani, Franosch 
and Frey. An extended model of the PFF model is used to describe an intra-cellular 
transport of the single-headed kinesin (KIFIA) motor [11]. 
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In the open boundary ASEP, we can draw the phase diagram by the incoming rate 
and outgoing rate at the boundaries. On the boundary between the low density phase 
and high density phase, it is known that the kink between a sparse region and a jammed 
region obeys Brownian motion. This kink motion in one-lane ASEP is also studied in 
terms of the domain wall theory [12, 13, 14] and the second class particle [15, 16, 17]. 
The kink motion in a two- lane TASEP is studied in ref.[18], and the kink motion in the 
PFF model is studied in ref.[19, 20, 21, 22]. It is notable that the kink is trapped by a 
harmonic potential in the PFF model [19, 20]. 

In this paper, we discuss the kink motion of TASEP with detachment and 
attachment. Our method can be used for any number of sites of detachment and the 
attachment. The effect of detachment of particles at the middle site of the system 
was studied in ref.[23] by dividing the system into two systems of TASEP, though the 
kink motion was not discussed there. We demonstrate that one site of detachment and 
attachment attracts the kink where attractive potential is a linear function of distance 
from the site of detachment and attachment. Even when we generalize the model with 
many sites of detachment and attachment, the kink feels the linear combination of the 
linear potentials for one-site of detachment and attachment. Though the analysis of 
linear combination of linear potential we derive the harmonic potential in PFF model. 

The organization of this paper is as follows. We first introduce the model with 
detachment and attachment at one site and briefly review previous studies on the kink 
motion in TASEP without detachment and attachment in section 2. Wc present a theory 
of a kink motion in TASEP with detachment and attachment at one site, and extend 
our analysis to the kink motion in TASEP with detachment and attachment at many 
sites including the PFF model. In section 3, we compare our analysis with the Monte- 
Carlo simulation and confirm the quantitative agreement between the theory and the 
simulation. Finally, we give concluding remarks. 

2. The kink motion 

2.1. A model for one site of detachment and attachment 

Let us explain TASEP model with detachment and attachment of particles at one site 
under the open boundary conditions. TASEP is defined on a one-dimensional lattice of 
L sites. Each particle hops forward when the site in front is empty. Here we choose 
right as the drift direction of particles. The open boundary condition is specified by a 
particle attachment rate a at the left end of the system and a particle detachment rate 
P at the right end. Here, the particle is also detached by the rate and attached by 
the rate Wa at the site xq. TASEP with detachment at one site has first been introduced 
in the ref.[23]. They divide the system into two systems of TASEP by the site of the 
detachment, and introduce the effective hopping rate to connect the two systems at the 
site Xq which is fixed to the middle of the system to include the effect of detachment. 
However, we can not explain the kink localization near xq by introducing the effective 
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hopping rate between two systems because the kink moves under the reflective boundary 
condition of a virtual boundary in one subsystem where the kink cannot cross through 
the virtual boundary. Thus we need another method to explain the effect of detachment 
and attachment to the kink motion. 

The Brownian motion of the kink is characterized by two parameters, the drift 
velocity Vr and the diffusion constant Dt- In the case of TASEP, Vr and Dt are 
written as [17] 

Vr-1-X,-Xr , Dr- 2(A. - A,) ' 

where Xe is the density in the left of the kink and A^ is the density in the right of the 
kink which satisfy A^ < A^. Thus once A^ and A,- are known, the motion of the kink can 
be determined. 

It is known that the Brownian motion in a potential U (x) is written by the Langevin 
equation 

where ^(t) is the random force satisfies 

(cw) = o {at)cm-^mt): (3) 

where S{t) is Dirac's delta function. The corresponding Fokker-Plank equation to eq.(2) 
is 

dP(x,t) d fdUix) \ , 

^ - ' ^ ^ +D—]P{x,t), (4) 



Pstix) = Cexp 



(5) 



dt dx \ dx dx ^ 

where P{x^ t) is the probability distribution of the Brownian particle at position x at 
time t. Equation (4) has the steady solution Pst{x) 

U{x 

where C is the normalization constant. We expect that this well-known result can be 
used to describe the steady distribution of the kink in our model. 

For TASEP without detachment and attachment, A^ and A,, satisfy 

A^ = q; , A^ = 1 — a (6) 

ior a — P < 1/2. Prom eq.(l), Vr and Dt are respectively given by 

a(l — a) 

Vt^O , Dt= ; , ^ (7) 
i — Za 

If there is a site with detachment and attachment, the phase diagram may be 
modified from that of the original TASEP. However the difference is expected to be 
small for small rate of the detachment and attachment {wa ~ <^ 1). Thus the kink 
picture may be still vahd for a = /3 -C 1/2. 

Let us consider the kink motion in TASEP with detachment and attachment at 
one site. If there is no detachment and attachment of particles, A^ and A^ are given by 
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Figure 1. A schematic picture of the kink motion in TASEP with one site of 
detachment and attacliment.The figure (a) shows the case when xq is in the low density 
region of the kink, and the figure (b) shows the case when xq is in the high density 
region of the kink. The fine shows the density profile. 



eq.(6). However, as a result of detachment and attachment the density is deviated from 
eq.(6). Let us consider the deviation of the density for two cases :(a) the site of the 
detachment and attachment is in the low density region of the kink and when (b) the 
site of detachment and attachment is in the high density region of the kink. 

Let us think of the case (a). The density in the left of the kink position xq is not 
affected by detachment and attachment at xq but the density in the right of xq becomes 
higher as in the Fig. 1(a). This situation takes place because attachment process is likely 
for low-density region in the one-way motion of the particles. We introduce p+ as the 
density in the right of Xq when xq is in the low density region of the kink. Thus 
should be = p+ while A^ is unaffected as A,. = 1 — a when xq is in the low density 
region of the kink. 

When Xq is in the high density region of the kink as in the case (b), the density 
is affected as in the Fig. 1(b). In this case, density in the right of xq is not affected by 
detachment and attachment but the density in the left of Xq is changed by detachment 
and attachment. The high density of the particles induce the excess of the detachment 
current to the attachment current for Wa ~ Wcl- Thus A^ and A^ should be replaced by 
X(, = a and A,. = p_ when xq is in the high density region of the kink. 

Now, let us determine the value of p_ and p+. Since the exact solution in the steady 
state of this model is not known, we adopt the decouple approximation in which the 
mean currents Ji in the low density region of the kink and in the high density region 
of the kink are given by = A£(l — A^) and = Ar(l — A^) respectively. To estimate 
p+ and p_, we use the current conservation at the position xq. For the case when xq is 
in the low density region of the kink, the current from the left oi Xq is a{l — a) and the 
current to the right of Xq is p+(l — p+). The detach current from xq is WdP+ and the 
attach current at xq is Wa{^ ~ P+)- Thus the equation of the current balance at xq is 
given by 

a{l -a)+ Wa{l - p+) = WdP+ + p+(l - p+). (8) 
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1 + u'a + "V - J (1 + iv„ + »',/)- - l»v, - in (1 - a) 
P-f = ^ • (9) 

Here we use the condition that p+ should be reduced to (6) for Wa — wa — 0. Similarly, 
the current balance equation of p_ is written as 

- p_) + Wa{l - P-) = WdP- + ol{1 - a), (10) 

and its solution is given by 

^ _ 1 - - Wrf + ^/(l - Wg - WdY + 4wa - Aa{l - a) ^^^^ 

Prom eq.(l) with the consideration around Fig.l, the drift velocities V± are given 

by, 

y+ = a — p+ and V^ — l — a — p^. (12) 
Similarly the diffusion constants D± are given by 

+ 2(1-q;-p+) 2(p_-q;) ^ ' 

Here the quantities with the suffix +/— represent those in the right/left of Xq. Using 
\^ — a,\r — p- ior X < Xq and = p+, = 1 — a; for x > Xq, the steady solution of 
the Fokker-Planck equation can be written as, 



P{x) = C'exp 



— (a; - xq)6{xq - x) + 7;— (a; - xq)6{x - xq) 



(14) 



where C is a normalization constant, and ^(x) = 1 for x > and 9{x) — otherwise. 
For w — Wa — Wd, the relations 

p-^l-p+ , \V+\^\V_\ , D+^D_ (15) 

are satisfied. Thus the solution (14) becomes 

P{x) = C"e-^l^-^°l, (16) 

where V = \V-\ = \V-^\ and D = = From the comparison of eq.(5) with eq.(16), 
the potential energy U (x) is given by 

U{x) = V\x~xo\. (17) 

This is one of the main results in this paper. The validity of our theory will be confirmed 
by Monte- Carlo simulation in section 3. 

2.2. Many sites of detachment and attachment 

Here we generalize the model in the previous subsection to a model with many sites 
of detachment and attachment. To ensure the domain wall picture, we assume that 
w = Wa = Wd are small and the relations eq.(15) are satisfied. Substituting the expansion 
p_i_ by w that 

p+ ~ q; + w H , (18) 
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we obtain 

V^w (19) 

In this linear regime, we can obtain the probabihty function of the kink in many sites 
of detachment and attachment. 

Let us consider the system with N sites of detachment and attachment. Thus the 
system is divided hy N + 1 segments in which each segment is bounded by the sites of 
detachment and attachment or the boundaries. 

When the kink is located in the jth segment, the number of sites of detachment and 
attachment in the left of the kink is j — 1. For each site of detachment and attachment, 
the current should be conserved. The density increases by w in the low density region 
as in the eq.(18) in each segment. Thus the density Xe in the j-th segment becomes 
Xi = a+w{j — 1). On the other hand, the number of sites of detachment and attachment 
in the right of the kink is N — j + 1. Thus the density in the j-th segment is 
Xr — 1 — a — w{N — j + 1). Thus, the drift velocity Vj and the diffusion constant Dj in 
the segment j are respectively given by 

Vj^w{N-2{j-l)}, (20) 

and 

^ 2a{l - g) + (1 - 2a)wN + w\{N - j f + f) 
^ 2{l-2a-wN) ' ^ ' 

Let us omit j-dependence in the eq.(21) by neglecting the term proportional to w"^. 
Thus the diffusion coefficient is reduced to 

_ 2a(l-n) + (l-2a)».-iV 

2{l - 2a -wN) ' ^^^^ 

which is independent of j. Prom eq.(5), the kink probability distribution Pj{x) in the 
segment j is given by 

2w{N - 2{j - !))(! -2a- wN) 



Pj (x) = Cj exp 



{X - Xj) 



(23) 



2a{l - a) + (1 - 2a)wN 
where Cj is a constant determined from the normalization and the compatibility relation 

Pj{xJ)^PJ+^{xJ). (24) 

Thus we obtain the formula applicable to any number of sites of detachment and 
attachment. It should be noted that the denominator in the exponential function in 
eq.(23) contain wN. This is because we simply put D in eq.(22) into eq.(5). It may be 
controversial to contain wN in the denominator as the systematic approximation, but 

the expression gives us the accurate result as will be shown later. 

This method is also applicable to the PFF model [10], where the number of 
detachment and attachment sites N is N = L, the segment length is 1, and the position 
of the site of detachment and attachment Xj is Xj = j. Thus the distribution (23) is 
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reduced to Pj{j) — Cj. We write Pj — Pj{j) for the simplification. The compatibihty 
relation (24) becomes 

-w{L^2{j-l)y 



Cj = Cj_i cxp 

and this recursion relation gives 

i-i 



D 



Cj = Ci exp 



m=l 



Thus the distribution is given by 



Pj — Ci exp 



w 
D 



(25) 



(26) 



(27) 



2 ' 4 

This result can be also derived by the superposition of the potential. The harmonic 
potential is realized by the superposition of the potential (17) as 

(28) 



UpffU) = w\j -xo\c^w (j - ^^Y^ ) + 

XO = l ^ ^ 



Thus the distribution of the kink position is given by 
PpffU) = C" exp 



W . . L + 1, 2 



(29) 



where D is given by eq.(22), and the probability distribution function of the kink position 
(29) is identical to eq.(27). 



3. Simulations 



Now let us check the quantitative accuracy of our theoretical argument. Wc compare 
our analysis with the results of Monte-Carlo simulations. The simulation is carried out 
by the random update scheme [3, 24], which is realized by choosing the bond between 
two neighbouring site randomly and move the particle in the chosen sites stochastically. 

We use the motion of the second class particle which is a tracer particle of the kink 
position to detect the kink position [15, 16, 17]. If we write for a hole (empty site), 
1 for a particle and 2 for the second class particle, the second class particle moves as 
(2,0) (0,2) and (1,2) (2, 1). The hole is moved to the left of the second class 
particle, and the particle is moved to the right of the second class particle. Thus the 
second class particle is positioned between the low density region of particles and the 
high density region of particles. Thus the second class particle can detect the kink 
position. We compare the probability function of the kink position obtained by the 
simulation with the distribution function (14). 

As shown in Fig. 2, we obtain good agreement between the simulation and the 
theoretical results in eq.(14). We plot the results of our simulation by x and the 
solution of the Fokker-Plank equation by the solid line. The horizontal axis is the kink 
position and the vertical axis is the probability distribution function of the kink position 
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Figure 2. Comparison of the distribution of the kink position between the simulation 

(x) and the solution of the Fokker-Plank equation (solid line). The horizontal axis is 
the kink position and the vertical axis is the distribution function of the kink position 
in the steady state plotted in the semi- log scale. The parameters used are a = /3 = 0.1 
and Wa = Wd = 0.01 for (a), and a = p = 0.1, = 0.02 and Wd = 0.01 for (b). We 
set the system length L = 200 and the position of detachment and attachment occur 
at xo = 100. 
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Figure 3. Comparison of the probability distribution of the kink position between 
the simulation and the solution of the Fokker-Planck equation (solid lines) in eq.(23) 
for the case which has 2(x), 3(+) and 4(*) sites of detachment and attachment. 
The horizontal axis is the kink position and the vertical axis is the probability 
function of the kink position in the steady state plotted in the semi-log scale. The 
parameters used are a = /3 = 0.1 and w = 0.01, and the positions of detachment and 
attachment are xi = 50,^2 — 150 for N = 2, Xi = 50i,{i = 1,2,3) for A'' = 3 and 
Xi = 40i, {i = 1, 2, 3, 4). We set the system length L = 200. 

in the steady state. In Fig. 2 (a), the parameter is set to be Wa = and the distribution 
function is symmetric around x — xq. In Fig. 2(b), the parameter is set to be Wa 7^ Wd 
and the distribution function is asymmetric around x — Xq. In both cases, the boundary 
parameters are set tohe a — (3 — 0.1 and the system length L is L = 200. In Fig. 3, we 
compare the theoretical results with the results of our simulation in cases of 2, 3 and 4 
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Figure 4. Comparison of the probability distribution of the kink position between 
the simulation (x) and the solution of (eq.(23))the Fokker-Planck equation (solid line) 
of the PFF model. The horizontal axis is the kink position and the vertical axis is the 
distribution function of the kink position in the steady state plotted in the semi-log 
scale. The parameters used are a = (3 = 0.1,w = 0.0001 and L = 200. 



sites of detachment and attachment. The parameters used in this case are a = (3 = 0.1 
and w = 0.01, the position of detachment and attachment are Xi = 50, 0:2 = 150 for 
N = 2, Xi = 50i, {i = 1,2,3) for = 3 and Xi = 40i, {i = 1,2,3,4). The results of 
our simulation are plotted by x for A?^ = 2, + for = 3 and * for = 4, and the 
theoretical results are plotted in the solid lines. In Fig.4, we compare the theoretical 
results with the simulation results in the PFF model. The parameters at the boundaries 
are a — P — 0.1 and w — 0.0001. The results of our simulation give quite well agreement 
of the theoretical prediction in eq.(23) in all cases. 

4. Conclusion 

We have demonstrated that the kink motion in TASEP with detachment and attachment 
can be described by the Brownian motion under the influence of the attractive force from 
detachment and attachment sites. We have obtained the attractive potential to the kink 
and the diffusion constant of the kink. We demonstrate that the superposition of the 
potentials of our model gives good results for any number of sites of detachment and 
attachment when the rates of detachment and attachment are small. We compare our 
result with the simulation and have confirmed that our theoretical prediction of our 
theory gives quantitatively correct results. We also explain the reason why the kink in 
PFF model feels a harmonic potential [19, 20], and succeed the quantitative estimation 
of the potential. 
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